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Coefficient of the x? term in the equation of a 
line on a constant life fatigue diagram where 
minimum,stress is x and maximum stress is jy. 

(R = Ax” + Bx + C - y) 

An assigned value of tl.or -l. 

A stress used in the calculation of plastic strain. 
The absolute value of DIF. 

The absolute value of ASMAX or of ASMIN, as assigned. 
The absolute value of ASMAX. 

The absolute value of ASMEAN 

The absolute value of ASMIN. 

The absolute value of R(4). 

The absolute value of R(7). 


The name of a routine calling for the absolute value 
of a quantity. 


Stress concentration factor, K,- 
The product (AKT) (STMAX). 

The quantity (ASMAX + ASMIN)/2. 

The product (AKT) (STMIN). 

Average value of SIGMIN over an interval. 

Average value of SIGMAX over an interval. 
Coefficient of the x term. (see A) 

A stress used in the calculation of plastic strain. 


The constant. (see A) 


Inverse of the Coffin-Manson slope. 


The number of cycles in an interval. 


CYCLES 


ENEP 


ENN 
ENNCYC 


EPSD 
EQRES 
EX 


EXP 


EXPO 


El 


The calculated number of cycles expected to be 
indicated on a constant life fatigue diagram for the 
applied combination of maximum and minimum stress. 
The Residual Stress Relaxation Constant. (see ENEP) 
Damage. 

Dev:imal or real value of integer K after conversion. 


A portion of a least-squares-method solution. 


The difference between residual stress and equilibrium 
residual stress. (RES(I) - EQRES) 


A portion of a least-squares-method solution. 

A variable used in the calculation of the number of 
cycles to be considered as an interval for relaxa- 
tion determination. 


The elastic modulus. 


The number of cycles from the beginning of the relax- 
ation process to the end of the current interval. 


The number of cycles required for overload residual 
stress effect to return to within one-tenth of its 
original difference from equilibrium conditions. 
(Nop * c1/(aBM)™! (apmean) ©?) 

The number of applied cycles at a load level. 


The ratio of the number of applied cycles to the 
number of cycles to failure. (ENN/CYCLES) 


LCF strain intercept. 
Equilibrium residual stress. 


An exponential function depicting the relaxation 
of residual stress. 


The name of a routine calling for the exponential 
value of a quantity. 


An exponent. The power of 10 which indicates the 
number of cycles to failure. 


Residual stress relaxation exponent. 


Residual stress relaxation exponent. 
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FLOAT 


I 
IBLOCK 


IFIX 


IN 


IPRINT 
IRAIN 
IRPCM 


1STEP 


ITYPE 


The name of a routine calling for integer-to-real 
conversion. 


A variable subscript. 


The identifying number of a block, the blocks being 
numbered consecutively from 1 to NBLOCK. 


The name of a routine calling for real-to-integer 
conversion. 


The number of steps input to the range-pair counting 
subroutine. 


Value controlling the WRITE statements. 
A counter. 


Value controlling entry into the range-pair counting 
subroutine. 


The identifying step number, the steps being numbered 
from 1 to NLEVEL. 


The identifying type number, the types being numbered 
from 1 to NTYPE. 


A variable subscript. 

Value of +1 or 0, as assigned for branch determination. 
Value of -1 or 0, as assigned for branch determination. 
An index variable. 

An index variable. 

An index variable. 

An index variable. 

An index variable. 


The number of steps output from the range-pair 
counting subroutine. 


An index variable. 
An index variable. 
An index variable. 
An index variable with values of N=4-7 indicating 


the power of 10, and thus identifying a particular 
life cycle curve. 
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NBLOCK 


NDECK 
NFLAG 
NFLAG2 
NLEVEL 


RES 
RNCYC 


SIGMAX 
SIGMIN 
STMAX 
STMIN 
SUMDEL 
SUMENN 


SUMNC 
SUMR 
SUMRN 
SUMR2 
SUMR2N 
SUMR3 
SUMR4 
TITLE1 


The total number of times to execute a block of 
loads. 


The number of data decks to be run sequentially. 
An integer used as a counter. 
An integer used as a counter. 


The total number of steps, or levels, of loads in 
a block. 


A subscripted variable used to indicate which types 
of loads are experienced in which blocks. 


The total number of different types. 
Plastic strain. 
Residue term in damage calculation. 
Residual stress. 


The number of cycles for a level after exitting 
the range-pair counting subroutine. 


Maximum stress. 

Minimum stress. 

Maximum applied stress. 

Minimum applied stress. 

Summation of damage for a flight. 


Accumulated total of applied cycles. (Summation of 
ENN) . 


Accumulated cycle ratio. (Summation of ENN/CYCLES) 


Summation of R(N), N=4,7. 

Summation of NR(N), N=4,7. 
Summation of R(N)*, N#4,7. 
Summation of NR(N)2, N=4,7. 
Summation of R(N)°, N=4,7. 
Summation of R(N)*, N=4,7. 


Identification of the source of the S-N data. 


TLE Tensile limit load. 


T1,TM2 Material type. 
TTYS One-fifth of tensile yield stress. 
Ts Tensile yield stress. 


TL, ¥2 Test identifying information. 


X Variable equivalent to SIGMIN. 
x Variable equivalent to SIGMAX. 
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I. INTRODUCTION 


In an aircraft fatigue life monitoring program, deter- 
mination of the actual loading environment which produces 
the fatigue is a major problem. With both military trainer 
and fighter aircraft, which may be used for a variety of 
duties, regular or continuous load recording programs have 
to be considered mandatory. Such load recording programs 
provide the data to calculate the consumed fatigue life of 
operating aircraft and are used in selecting design spectra 
for future aircraft designs. 

To facilitate the calculations of fatigue life from 
such data, it is important to know not only the loading 
magnitudes but also the time sequence of the loading. The 
fatigue life of a local critical point varies due to 
residual stresses remaining after the application of a load 
causing local plasticity in tension or compression. For 
instance, a peak tensile load into the plastic flow region 
leaves a residual compressive stress, which lowers the 
local magnitude of the following tensile stresses, thus 
increasing the fatigue life. In a similar manner a peak 
compressive load into the plastic flow region, if it were 
to occur, would leave a residual tensile stress, which 
decreases the fatigue life. 

Currently, fatigue monitoring of naval aircraft is 


based on the total number of g readings recorded at four 
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selected levels by an exceedence level counting accelerometer. 
Using micropressors in the near future it will be possible 

to record each maximum and minimum load level experienced 

by an aircraft in sequence. The data collected can be used 
to monitor the fatigue life of a structure via the deter- 
mination of damage accumulated at a point found to be 
critical in a structural test of a prototype. It should 

be noted that this calculation uses a theoretical model. 
Damage is not observable or measurable. 

The objective of this thesis is to use a computer program, 
employing the range-pair counting technique on time sequence 
recorded maximum and minimum aircraft loads and a relaxation 
model to consider residual stresses, to calculate plastic, 


elastic, and total aircraft structural damage. 


II. GENERAL DESCRIPTION OF THE COMPUTER PROGRAM 


The computer program used in this fatigue analysis is 
divided into four modules. Each module is clearly labeled 
in the program listing (page 93). In module I the input 
data are read and assimilated in preparation for future 
calculations. The input parameter requirements are presented 
in detail beginning on page 18, and six sample sets of 
input data are illustrated starting on page 40. 

When the load sequence randomization technique is used 
to input loading data, it is a part of module I. The load 
sequence randomization technique uses the computer library 
subroutine RANDU to place 10 percent of the MIL spectrum 
A positive loads (Table I) in a random order. Each of the 
positive loads is paired with a minimum load of 11 percent 
limit load, or l-g. During the randomization each load has 
an equal probability of selection. A counter restricts the 
number of times a value is selected to the number of 
occurrences of the particular load level in MIL spectrum A. 

Two variations on the randomization of the MIL spectrum 
A load levels are available. The first option permits 
randomization of the negative MIL spectrum A load levels 
using the same technique as in the case of the positive 
MIL spectrum A load levels. Each negative load level is 
randomly paired with a positive load level. Since there 


are not as many negative load levels in MIL spectrum A as 


TABLE I 


FREQUENCY OF MANEUVER LOADS 


NUMBER OF TIMES PER THOUSAND HOURS THAT 
LOAD FACTOR IS EXPERIENCED 


PERCENT OF MAXIMUM (POSITIVE) FLIGHT MANEUVER LOAD 
SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 

35 17000 

45 9500 

55 6500 

65 4500 

75 2500 

85 1360 

95 440 

105 150 

115 40 

125 16 

TOTAL ~32006 
PERCENT OF MAXIMUM (NEGATIVE) FLIGHT MANEUVER LOAD 

SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 

0 500 

10 200 

20 100 

30 60 

40 35 

50 30 

60 25 

70 20 

80 15 

90 10 

100 5 

110 3 

TOTAL 003 
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there are positive load levels, the excess load levels are 
paired with 1l-g loads. Mixing of the negative MIL spectrum 
A load levels and the 1-g load levels is accomplished using 
the computer library subroutine RANDU. For this ordering 
thereis an 80 percent chance that each negative load selected 
will be a 1l-g load. A counter ensures the proper numbers of 
each load type are included. Eighty percent probability of 
selection of a 1-g load was used to spread the smaller 
number of negative MIL spectrum A loads throughout the 
sequence, and also to speed computer operation by not having 
the matrix of negative MIL spectrum A loads addressed so 
often, after the counter indicated the matrix elements had 
all been used. 

The second option for randomizing loads provides for the 
inclusion of a ground cycle between each flight. The ground 
cycle loading is taken as -2500 psi, or -8 percent limit 
load. Eight percent is the level of the ground cycle load 
used by the A7 manufacturer. There are 4201 positive 
loading events in 100 hours of MIL spectrum A flight time. 
Each flight is arbitrarily taken as one hour in length, so 
a ground load is paired with every forty-second positive 
loading event. 

In module II the local stresses and strains are deter- 


mined from the input data. If the local peak stress is in 


: 
% 
1 
i 
4 
: 
5 
g 
= 
> 
i 
Fe 
z 
& 


the plastic region, the specimen is assumed to unload 
elastically, leaving a local residual compressive stress 


in the material. In the analysis developed by Potter [Ref. 6], 
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, the transient portion of the residual stress so produced 
relaxes toward zero, or an equilibrium residual stress. 
After the stresses and strains are calculated, the local 
stress cycles are counted in module III using the range- 
pair counting technique. In module IV the fatigue damage 
is calculated. Damage is determined separately for elastic 
and plastic strain events. The elastic and plastic damage 
for each flight is listed as output, as well as the total 


cumulative damage from any previous flights. 


A. INPUT DATA REQUIREMENTS (MODULE I) 

Samples of data are listed starting on page 40; format- 
ting is described beginning on page 93. Both pages should 
be consulted while reading the following descriptions. 

1. Data Card 1 

Only one of these cards is required or permitted 
for each program run. 

NDECK = the number of data decks to be run sequen- 
tially. The input on this card will be common to all data 
decks run. It is not necessary for the different data 
decks to have any parameters in common. 

IPRINT = the value controlling the write statements. 
Output available during the process of the analysis 
includes: 

a. Maximum and minimum applied stress of each cycle 
and the local stress response throughout the spectrum. Also 
printed out is the residual stress, equilibrium stress, 


applied cycles, and the equilibrium period for each cycle. 
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b. The elastic local stress history as input into 
the range-pair subroutine and the resulting range-paired 
spectrum. 

c. The maximum plastic strain occurrence during the 
spectrum and the damage associated with each strain reversal. 

d. The accumulated damage associated with the plastic 
strains. 

e. The range-paired elastic stress spectrum and the 
damage associated with each level for block loading, or for 
each cycle for RANDU generated loadings. 

f. The accumulated damage associated with the cur- 
rent block of loading, including the plastic strain damage, 
and the total damage since the initiation of cycling. 

The value assigned to IPRINT controls which pieces 
of dataare included in the printout. 

If IPRINT = 1, all six items listed above are 
printed for each flight or block of loads. 

If IPRINT = 2, all items listed above except b, 
are- printed. 

If IPRINT = 3, only items d and f listed above are 
printed. 

An example of each output option is given beginning 
on page . When the number of loads is more than fifty, 
the output is quite voluminous, and a paper usage default 
may interrupt program operation. With large numbers of 


loads, IPRINT should be set equal to 3. 


IRPCM = the value controlling entry into the range- 
pair counting subroutine. 

If IRPCM = 1, the range-pair counting method is used. 

If IRPCM = 2, the range-pair counting method is 
skipped. 

Each data deck run sequentially must contain the 
following cards. The number of data decks was indicated by 
the value assigned NDECK. 

2. Gar i 

Information on this card is descriptive, alphanumeric 

data used only as a heading test identification. 
3. Card 2 

The name of the material type and four material con- 
stants are listed on this card. The material type is an 
alphanumeric entry. The tensile yield stress, LCF strain 
intercept, inverse of the Coffin-Manson slope, and the 
modulus of elasticity are real number entries. The LCF 
strain intercept is the numerical value of the ordinate 
intercept in the log-log Coffin-Manson plot of plastic 
strain range vs. cycles to failure. Units of tensile yield 
stress and modulus of elasticity are ksi. 

4. Cards 3 Through 6 

Each card contains the three coefficients of a second 
order least-squares curve fit for S-N data obtained from a 
Goodman diagram published in MIL-HDBK-S5A. The curves are 
fit for lives of 104, 10°, 10°, and 10’ cycles. Space is 


provided on each card for an alphanumeric indication of 
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the data source. Only the data source printed on card 6 


will be retained in computer storage for subsequent printout. 
Columns 72 through 80 are used for information only and are 
not read by the computer. 
S. Card 7 
Constants to be used in the equilibrium period cal- 
culations are listed on this card. 
6. Card 8 
The stress concentration factor is the only entry 
on this card. 
7. Card 9 
The first entry on this card is the number of blocks 
of loads, or the number of flights. The terms block and 
flight are used interchangeably. This entry is the number 
of times the list of loads is to be repeated. The next 
entry is the number of loads. The number of loads will 
always be the same as the quantity of cards 11 used as input. 
The last entry is the number of types of loads. Loads are 
indicated as being a different type in order to alter the 
loading pattern during subsequent flights. It is important 
to understand that the program will always consider load 


types in ascending numerical order. A flight listing load 


types in the order 1, 2, 3, 4 will produce identical results 
with a flight listing load types in the order, 4, 3, 2, 1. 
8. Card 10 
Only the limit load is listed here, with units of 
ksi. The maxima and minima of all cycles are input as a 
decimal fraction of this load entry. 
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9. Cards ti 

Information identifying each load cycle is presented 
on these cards. The first entry is an integer representing 
the step, or number of the load. This integer increments by 
one on each card, with the last card having a value equal to 
NLEVEL. Other entries on these cards are the type of load, 
the minimum and maximum decimal fraction of the limit load 
for the particular step, and the number of consecutive 
cycles at this load level. It is important to notice that 
the format requires the number of cycles to be entered as 
a real number, although fractions of cycles will not be 
processed. 

10. Cards NLEVEL + 11 

The first entry is an integer representing the number 
of the block, or flight. This entry is a dummy variable and 
is not used in any subsequent calculations. The other 
entries on this card indicate which types of loads are to 
be processed on this particular flight. Each type of load 
may be processed only once on each flight. The order in 
which the types are listed is of no significance, since the 
load-types are automatically processed in ascending numerical 


order. 


B. RANGE-PAIR COUNTING METHOD 

In assessing the life consumed by an individual aircraft, 
damage is calculated from Miner's Law using S-N data. This 
type of life calculation is of greater value if the cycle 


counting method used takes into account the actual load-time 


history of the aircraft. Many load cycle counting methods 
have been developed and used. Fatigue test experience indi- 
cates that a useful counting technique for aircraft fatigue 
must take into account the loading sequence as well as the 
magnitude of the maximum and minimum load peaks [Ref. 10]. 
The two most realistic counting techniques are the range- 
pair counting method and the rain-flow counting method. 

Both methods experimentally yield approximately the same 
results. The major difference between the two methods is 
that the rain-flow method counts in terms of load ranges, 

or half cycles, whereas the range-pair method records the 
time history in terms of complete load cycles. When counting 
in terms of full cycles, the two methods are equivalent. 

The merit of counting half cycles is not important when 
analyzing aircraft fatigue because the large number of load 
reversals does not permit a half cycle to influence the 
numerical results. The range-pair counting method corresponds 
to the stable cyclic stress-strain behavior of a material 

in that strain ranges counted as cycles will form closed 
stress-strain hysterisis loops [Ref. 4]. 

The computer program employed in this thesis has the 
option of using the range-pair counting method, or simply 
using each peak in the order that it occurs. When the range- 
pair counting method is not used, only momentary load values 
are considered, and information regarding the cyclic stress- 
strain pattern, which is important in fatigue calculations, 


is lost. Also, minor load variations, which are of little 
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importance in fatigue calculations, are counted as additional 


cycles. 

The range-pair counting method is illustrated in Figure 1. 
It counts a strain range as a cycle if it can be paired with 
a subsequent straining of equal magnitude in the opposite 
direction. For a complicated load history, some of the 
ranges counted as cycles will be simple ranges, such as 2-3, 
during which the strain does not change direction, but others, 
such as 1-8, will be interrupted by smaller ranges which will 
also be counted as cycles. In Figure 1 ranges are marked with 
solid lines and the paired ranges with dashed lines. 

Each peak is taken in order as the initial peak of a 
range, except that a peak is skipped if the part of the 
history immediately following it has already been paired with 
a previously counted range. If the initial peak of a range 
is a minimum, a cycle is counted between this minimum and 
the most positive maximum which occurs before the strain 
becomes more negative than the initial peak of the range. 

For example, in Figure 1 a cycle is counted between peak 1 
and peak 8, peak 8 being the most positive maximum that 
occurs before the strain becomes more negative than peak l. 
If the initial peak of a range is a maximum, a cycle is 
counted between this maximum and the most negative minimum 
which occurs before the strain becomes more positive than 
the initial peak of the range. For example, in Figure 1 a 
cycle is counted between peak 2 and peak 3, peak 3 being 


the most negative minimum before the strain becomes more 


Strain —~ 


Figure 1 


Illustration of Range-Pair Counting Method 
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positive than peak 2. Each range that is counted is paired 
with the next straining of equal magnitude in the opposite 
direction, explaining why complete cycle rather than half 
cycle counts are made. For example, in Figure 1 part of the 
range between peaks 8 and 9 is paired with the range counted 


between peaks 1 and 8. 


C. DAMAGE CALCULATION 


Damage due to plastic strain and damage due to elastic 
strain is calculated separately. The behavior of the specimen 
is assumed to be elastic-perfectly plastic as illustrated in 
Figure 2. Behavior in compression is assumed to be the same 
as the behavior in tension. The total damage is determined 
using the Miner's summation D.= Zn/N_ where 

D = total damage 
n = number of cycles at a particular load level 


N = number of cycles to failure at a particular 
load level 


Failure of the specimen is defined as occuring when D equals 
unity. 

Plastic strain damage is calculated using the Coffin- 
Manson theory [Ref. 5]. The Coffin-Manson theory is based 
on an experimentally determined log-log plot of plastic 
strain range vs. cycles to failure as shown in Figure 3. 

The slope of the line in Figure 3 is approximately -1/2 
for all metals, with the ordinate intercept making the 
behavior of each metal unique. The equation used for cal- 


culation of plastic strain damage is 
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stress, 6 


tensile yield 
stress, Sty 


Strain, e 


Figure 2. ELASTIC-PERFECTLY PLASTIC 
BEHAVIOR 


log plastic j 
strain range, slope approximately -4 


®p 


log cycles to failure, N 


Figure 3. COFFIN-MANSON DIAGRAM 


27 


=1 


F=f ace" 


where N = number of cycles to failure at a specific plastic 


strain 
ep = magnitude of the plastic strain 
k = slope of the Coffin-Manson curve, approximately 


= 1/2 
c = ordinate intercept; experimentally for each metal 
type 
For strain levels in the plastic range, failure occurs in 
less than 104 cycles, and damage is calculated using the 
Coffin-Manson theory. For strain levels in the elastic range, 
the Goodman diagram is the basis for calculation of damage. 
The maximum and minimum local strain levels are sequentially 
compared with second order least-squares curve fit data from 
the Goodman diagram. The curve fits used are for lives of 
5 6 7 
,’ 


10°, and 10° cycles [Ref. 3]. The elastic damage 


ig”, 38 
is totaled using the Miner summation and combined with the 
plastic damage, computed in a similar way, to yield the total 
damage. The amount of plastic damage and the total damage 
for a particular flight, and the total damage including all 


previous flights is printed at the end of each flight. 


ee 


III. DISCUSSION OF RESULTS 


The load spectrum used in this analysis is the load spectrum 
A from MIL-HDBK-SA, and the loads are assumed to be applied at 
wing station 32 on an A7 aircraft wing [Ref. 2]. Load data 
were input as block loads and also as single randomized cycles. 
The Miner's Law summation of fatigue damage for this loading 
spectrum is .1300 per 1000 flight hours. Results obtained 
using this computer program are given in Tables II, III, and 
IV. It is important to notice that when fatigue damage is 
calculated without using the range-pair counting method, it 
makes only a minor difference on the total fatigue damage when 
the sequence of loading is changed 01 when negative loads are 
included. 

Block loads were arranged in LO-HI, HI-LO, and HI-LO-HI 
sequences. The results indicate the LO-HI sequence to be 
more damaging than the HI-LO sequence, with damage due to 
the HI-LO-HI sequence falling between the other two block 
load sequences. Inthe HI-LO sequence the HI loads leave 
a local residual compressive stress in the material, and 
the LO loads do less damage. 
In the case of RANDU generated load sequences, identical 
starter integers were used in order to evaluate the effect 
on the damage of inclusion of MIL spectrum A negative loads 
and ground cycles. The more negative loads that are included, 
the worse is the fatigue damage. The increase in damage with 


more negative loads is caused by the negative loads decreasing 
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TABLE II 


FATIGUE DAMAGE DUE TO POSITIVE MILSPEC A 
BLOCK LOADING, NORMALIZED TO 1000 HOURS 


SEQUENCE OF 


BLOCK LOADS USING RPCM* NOT USING RPCM* 


LO-HI - 16780102 - 16780108 
HI-LO -069243014 -069243014 
HI-LO-HI - 10864216 - 10839599 


*® RANGE-PAIR COUNTING METHOD 


TABLE III 


FATIGUE DAMAGE DUE TO MILSPEC A 
RANDOMIZED LOAD INPUTS, USING RPCM, NORMALIZED TO 1000 HOURS 


INTEGERS USED 


TO START RANDU MINIMUM LOAD VALUES 


11 PERCENT NEGATIVE GROUND 

LIMIT LOAD MILSPEC A CYCLES 
83745,54711,54487 -039226338 -046267733 -060848035 
13547 ,66549,7 -095124505 .085063167 10429338 
9,583,4777 043321513 -043640956 -036781987 
48621,3,491 -018672124 .034077277 -035881512 
73,559,1001 -031819174 -042006299 037130679 
357,833,1 -037665060 -049336702 .050973855 
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TABLE IV 


FATIGUE DAMAGE DUE TO MILSPEC A RANDOMIZED LOAD 
INPUTS, NOT USING RPCM, NORMALIZED TO 1000 HOURS 


INTEGERS USED 


TO START RANDU MINIMUM LOAD VALUES 
11 PERCENT NEGATIVE GROUND 
LIMIT LOAD MILSPEC A CYCLES 
83745,54711,54487 -040602274 -047398917 .049124956 
13547 ,66549,7 -040137805 -047727041 -049612150 
9,583,4777 -041204430 -048346408 -049540661 
48621,3,491 -041160211 -049230121 -050957575 
73,559,1001 -040792227 .047043338 - 049311966 
357,833,1 ; -041056350 -049220584 -050481893 
} 
1 | 
d 
4"! 31 
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the local residual compressive stress caused by positive loads 


into the plastic zone. This permits the local positve elastic 
stresses to be more damaging. 

Six plots are presented, beginning on page 33, to illus- 
trate the pattern of the individual load randomization process. 
Each plot portrays the first two flight hours of the aircraft. 
The range of values in the columns of Table III are caused by 
such variations in the load randomization pattern. 

The fatigue damage resulting from randomization of loads 
on a cycle basis, as shown in Table III, is significantly 
‘different from the damage due to the HI-LO-HI block loads 
illustrated in Table II. Vought Aeronautics Division con- 
ducted a fatigue study using randomized flights. Each flight 
was based on a mission profile specified by the Air Force. 
Fatigue damage resulting from the randomized flights was in 
close agreement with fatigue damage resulting from HI-LO-HI 


block loading. 
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IV. RECOMMENDATIONS FOR FUTURE RESEARCH 


Subsequent work in this area should be directed toward 
refining the input load sequence. The Naval Air Development 
Center collects accelerometer data from individual aircraft 
on a monthly basis. These load data could be used in this 
computer program to see if randomizing the load sequence over 
the life of the structure provides a different picture of the 
fatigue damage than randomizing by the month. Also, research 
is needed to determine the magnitude and sequence of the 
ground cycles. In this analysis all ground cycles were con- 
sidered of equal magnitude and equal spacing, which is not 
the case with a real airplane. A landing load spectrum should 
be used to be more accurate. 

The computer program itself could be improved by a better 
understanding and definition of the relaxation phenomena. 
Experiments to quantify the relaxation parameters will enhance 
the program. 

The Air Force is monitoring crack growth as the basis for 
their fatigue damage assessment. It would be productive to 
compare the total fatigue damage calculated by the two 


methods. 
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